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Abstract. Consider a non-standard numeration system like the one built over 
the Fibonacci sequence where nonncgative integers are represented by words 
over {0, 1} without two consecutive 1. Given a set X of integers such that 
the language of their greedy representations in this system is accepted by a 
finite automaton, we consider the problem of deciding whether or not A" is a 
finite union of arithmetic progressions. We obtain a decision procedure for this 
problem, under some hypothesis about the considered numeration system. In 
a second part, we obtain an analogous decision result for a particular class of 
abstract numeration systems built on an infinite regular language. 

1. Introduction 

Definition 1. A positional numeration system is given by a (strictly) increasing 
sequence U — (E/i)i>o of integers such that Uq = 1 and Cy := sup i>0 |~(7i + i/[/i] is 
finite. Let Ay = {0, . . . , Cy — 1}. The greedy U -representation of a positive integer 
is the unique finite word rep^ (n) — wi ■ ■ ■ wq over Ay satisfying 



n 



n = Wi Ui, wi 7^ and Wi Ui < Ut+\, Vt = 0, . . . , I. 

i=0 i=0 

We set rep;y(0) to be the empty word e. A set X C N of integers is U -recognizable 
if the language repu(X) over Ay is regular (i.e., accepted by a deterministic finite 
automaton, DFA). If x = X£ ■ ■ ■ xq is a word over a finite alphabet of integers, then 
the U -numerical value of x is 

t 

v&\u(x) = y^^Xj Uj. 

i=Q 

Remark 2. Let x, y be two words over Ajj. As a consequence of the greediness of 
the representation, if xy is a greedy [/-representation and if the leftmost letter of 
y is not 0, then y is also a greedy [/-representation. Notice that for m,n € N, we 
have to < n if and only if rep^m) < gen r epj/(n) where < gen is the genealogical 
ordering over A v : words are ordered by increasing length and for words of same 
length, one uses the lexicographical ordering induced by the natural ordering of the 
digits in the alphabet Ajj. Recall that for two words x,y £ A* v of same length, x 
is lexicographically smaller than y if there exist w,x',y' £ A v and a, b £ Ajj such 
that x = wax' , y = wby' and a < b. 

For a positional numeration system U, it is natural to expect that N is U- 
recognizable. A necessary condition is that the sequence U satisfies a linear recur- 
rence relation [55] , 

Definition 3. A positional numeration system U = ([/s)i>o is said to be linear, 
if the sequence U satisfies a homogenous linear recurrence relation with integer 
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coefficients. Otherwise stated, there exist k > 1 and constant coefficients a±, . . . , a k 
such that for all i > 0, we have 

(1) U i+ k = aiUi+k-i H h a k Ui, with a±, . . . , a k G Z, a fc 7^ 0. 

We say that fc is the order of the recurrence relation. 

Example 4. Consider the sequence defined by Fo = 1, F\ = 2 and for all i > 0, 
-Fi+2 = -Fi+i + Fi. The Fibonacci (linear numeration) system is given by _F — 
(Fi)i>Q = (1,2,3,5,8,13,...). For instance, rep F (15) = 100010 and val F (101001) = 
13 + 5 + 1 = 19. 

In this paper, we mainly address the following decidability question and its 
extension to abstract numeration systems. 

Problem 1. Given a linear numeration system U and a set X C N such that 
repy(N) and repj/(X) are both recognized by (deterministic) finite automata. Is it 
decidable whether or not X is ultimately periodic, i.e., whether or not X is a finite 
union of arithmetic progressions ? 

Notice that the regularity of rep[/(N) ensures that there exists a set ICN such 
that repjj(X) is regular, see also Remark 127)1 

Ultimately periodic sets of integers play a special role. On the one hand such 
infinite sets are coded by a finite amount of information. On the other hand the 
celebrated Cobham's theorem asserts that these sets are the only sets that are 
recognizable in all integer base systems (T2]. This is the reason why they are also 
referred in the literature as recognizable sets of integers (the recognizability being 
in that case independent of the base). Moreover, Cobham's theorem has been 
extended to various situations and in particular, to numeration systems given by 
substitutions [T3] . 

If we restrict ourselves to the usual integer base b > 2 numeration system, 
referred to in the sequel as b-ary system, defined by U% — bUi-i for i > 1 and 
Uq = 1, several results are known. J. Honkala showed in that Problem [TJ turns 
out to be decidable. Let us also mention [T], where the number of states of the 
minimal automaton accepting numbers written in base b and divisible by d is given 
explicitly. J. -P. Allouche and J. Shallit ask in [2] if one can obtain a polynomial time 
decision procedure for 6-ary systems. Using the logic formalism of the Presburger 
arithmetic, a positive answer to this question is given by J. Leroux in [3D] even 
when considering subsets of Z d , d > 1. In dimension one, ultimately periodic sets 
are exactly the sets definable in the Presburger arithmetic (N, +). 

Note that A. Muchnik has shown that Problem[T]turns out to be decidable for any 
linear numeration system U for which both repjj(N) and addition are recognizable 
by automata [32j . But it is a difficult question to characterize numeration systems 
U for which addition is computable by finite automaton, i.e., for which the language 
/0™-\™Pu(x)\ Tepu ( x )\ 
J \ 0™-\™Pu(y)\ Tepu ly)\ I x,y,zeN, x + y = z, m= max |rep y (t)|l 

t\ rn-\re Pu (z)\ Tepu ( z) J te{,,v,z} j 

where leading zeroes are prefixed for padding the shorter components to obtain 
three words of the same length, is regular (see for instance in which the char- 

acteristic polynomial of the sequence U is the minimal polynomial of a Pisot num- 
ber). In [50], the sequentiality of the successor function (i.e., the action of adding 1) 
is studied. If addition is computable by a finite automaton, so is the successor func- 
tion, but the converse does not hold in general. In particular, some examples of lin- 
ear numeration systems for which addition is not computable by a finite automaton 
are given in [5D]: for instance, the sequence defined by Ui = 3Ui-i + 2Ut-2 + 3Z7i_4 
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with integer initial conditions 1 = Uq < U\ < Ui < U3. So the decision techniques 
from 30, 32 cannot be applied to that system. Nevertheless, as we will see in Ex- 
ample [32l our decision procedure can be applied to this system. Notice also that in 
the extended framework of abstract numeration systems, one can exhibit systems 
such that multiplication by a constant does not preserve recognizability [HI [55] ISS] . 
Therefore the powerful tools from logic discussed above cannot be applied in that 
context either. 

The question studied in this paper was raised by J. Sakarovitch during the 
"Journees de Numeration!'' in Graz, May 2007. The question was initially asked 
for a larger class of systems than the one treated here, namely for any abstract 
numeration systems defined on an infinite regular language [29]. A shorter and 
partial version of this paper has been presented to the 33rd MFCS symposium in 
Toruh, August 2008 0. 

Observe that this decision problem for all abstract numeration systems is equiv- 
alent to the famous HD0L periodicity problem: given a morphism / and a coding 
g, decide if the infinite word g(/"(a)) is ultimately periodic, see [23 EE]- For the 
restricted case of the D0L periodicity problem, where only the morphism / is con- 
sidered, decision procedures are well-known 2 1 33 . Finally, questions connected 
to the ones addressed here have independently and recently gained interest [5] . In 
particular, a simple proof of Honkala's original result based on the construction 
of some automata is given in that paper. As for the logical approach considered 
by Muchnik and Leroux, the arguments given in [3] rely on the recognizability of 
addition by automata (which can be done for the classical fc-ary system but not 
necessarily for an arbitrary linear numeration system). 

The structure of this paper is the same as [26) . First we give an upper bound 
on the admissible periods of a [/-recognizable set X when it is assumed to be ul- 
timately periodic. Then an upper bound on the admissible preperiods is obtained. 
These bounds depend essentially on the number of states of the (minimal) automa- 
ton recognizing repu(X). Finally, finitely many such periods and preperiods have 
to be checked. For each of them, we have to build an automaton accepting the 
corresponding ultimately periodic sets (this implies that N has to be recognizable). 

Though the structure is the same, our arguments and techniques are quite differ- 
ent from [26] , They rely on the study of the quantity Njj{m) defined as the number 
of residue classes that appear infinitely often in the sequence (Ui mod m)i>o- Our 
main result — Theorem [22] — can be stated as follows. 

Theorem. Let U = (Ui)i>o be a linear numeration system such that N is U- 
recognizable and satisfying liirii_ ++00 f7i+i — Ui = +00. If lim m _^ +00 Nu(m) = +00, 
then Problem [1] is decidable. 

Actually our techniques cannot be applied to 6-ary systems (see Remark |2~5]) , 
which is the case treated by Honkala [25] , because in that case Nu(m) -f* +00 as 
m — ► +00. 

In Section [21 we give a characterization of the linear numeration systems U such 
that lim m ^ +00 Nu{m) = +00. To do so we use p-adic methods leading to a study 
of the sequence (Ui mod p v )i>o for all v > 1, where p is a prime dividing ak- 
in the last section, we consider again the same decision problem but restated 
in the framework of abstract numeration systems [29j . The definition of these 
systems is given in Section [H We apply successfully the same kind of techniques to 
a large class of abstract numeration systems (for instance, an example consisting 
of two copies of the Fibonacci system is considered). The corresponding decision 
procedure is given by Theorem 1421 As explained above, this result provides a 
decision procedure for specific instances of the HD0L periodicity problem. 
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All along the paper, we try whenever it is possible to state results in their most 
general form, even if later on we have to restrict ourselves to particular cases. 
For instance, results about the admissible preperiods do not require any particular 
assumption on the numeration system except linearity. 

2. A Decision Procedure for a Class of Linear Numeration Systems 

We will often consider positional numeration systems U = (Ui)i>o satisfying the 
following condition: 

(2) lim U i+1 -Ui = +oo. 

i — > + OG 

Notice that it is a weak requirement. Usually, the sequence U has an exponential 
growth, Ui ~ (3 l for some f3 > 1, and therefore ^ is trivially satisfied. (It is for 
instance the case for the numeration systems considered in Remark 1271 ) 

Lemma 5. Let U = (Ui)i>o be a positional numeration system satisfying ((2]). Then 
for all j , there exists L such that for all I > L, 

1Q i- 1 re P[ , (t) | vepu ^ £ = 0, . . . , JJj — 1 

are greedy U -representations. Otherwise stated, if w is a greedy U -representation, 
then for r large enough, I0 r w is also a greedy U -representation. 

Proof. Notice that rep u (Uj — 1) is the greatest word of length j in rep [/ (N), since 
rep [/ (C/j) = 1CP. By hypothesis, there exists L such that for all £ > L, Ue+i — Ui> 
Uj - 1. Therefore, for all I > L, 

lO^rep^-1) 

is the greedy [/-representation of Ut + Uj — 1 < Ut+\ and the conclusion follows. □ 

Example 6. Consider the positional numeration system Uo = 1, U\ = 2, U% = 3 
and Uzi+r = 3 i+1 + r for all i > 1 and r 6 {0, 1, 2}. This system does not satisfy 
dH) because Ui+i — Ui = 1 for infinitely many i's. We have rep [/ (2) = 10, but 
one can notice that for i > 0, 10 3z+1 10 is not a greedy [/-representation. Indeed, 
val[/(10 3i+1 10) = [/ 3 ( i+ i) + 2 = U 3i+5 has 10 3i+5 as greedy [/-representation. 

Remark 7. In the above lemma, one cannot exchange the order of the quantifiers 
about j and L. As an example, consider the sequence Ui = (i+l)(i+2)/2 for all i > 
0. This sequence satisfies the linear recurrence relation [/i+3 = 3[/i+2 — 3[/i+i + Ui 
and also Observe that, for all i > 1, k = Ui — 1 is the unique value such that 
Ui = Uk — C/fc-i because for all j > 1, Uj — Uj-i = j + 1. For all i > 1, 10 J is a 
greedy [/-representation and the greedy [/-representations of the form 10™10 i are 
exactly those for which n > Ui — i — 1. 

Remark 8. Bertrand numeration systems associated with a real number [3 > 1 
are defined as follows. Let Ap — {0, . . . , \(3~\ — 1}. Any x S [0, 1] can be written as 

x = Ci , with Ci € Ap 
i=i 

and the sequence (cj)i>i is said to be a (3 -representation of x. The maximal (3- 
representation of x for the lexicographical order is denoted dp(x) and is called the 
(3 -development of x (for details see [SU Chap. 8]). We say that a /3-development 
(ci)i>i is finite if there exists N such that Cj = for all i > N. If there exists m > 1 
such that dp(l) = t\'--t m with t m ^ 0, we set d%(l) :— (t% ■ ■ -t m ^i(t m — l)) w , 
otherwise ^(1) is infinite and we set dj 3 (l) := dp(l). 
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We can now define a positional numeration system Up = ([/i)i>o associated with 
/3 (see H). If dj(l) = {ti)i>i, then 

(3) U = l and Vi > 1, Ui = t 1 U l -i H h UU Q + 1. 

If /3 is a Parry number (i.e., dp(l) is finite or ultimately periodic) then one can 
derive from §3§ that the sequence Up satisfies a linear recurrence relation and as 
a consequence of Bertrand's theorem (lj linking greedy [/^-representations and fi- 
nite factors occurring in /3-developments, the language rep^ (N) of the greedy Up- 
representations is regular. The automaton accepting these representations is well- 
known |22j and has a special form (all states — except for a sink — are final and 
from all these states, an edge of label goes back to the initial state). We therefore 
have the following property which is much stronger than the previous lemma. If x 
and y are greedy [/^-representations then xOy is also a greedy [/^-representation. 

Example 9. The Fibonacci system is the Bertrand system associated with the 
golden ratio (1 + V5)/2. Since greedy representations in the Fibonacci system are 
the words not containing two consecutive l's [41] . then for i,t/£ rep F (N), we have 
xOy e rep F (N). 

Definition 10. Let X C N be a set of integers. The characteristic word of X is 
an infinite word x§X\Xi ■ ■ ■ over {0, 1} defined by Xi = 1 if and only if i 6 X. 

Consider for now X C N to be an ultimately periodic set. The characteristic 
word of X is therefore an infinite word over {0,1} of the form 

XqX\X2 ■ ■ ■ — fin" 

where u and v are chosen of minimal length. We say that the length |tt| of u (resp. 
the length \v\ of v) is the preperiod (resp. period) of X. Hence, for all n > 
n e X if and only if n + \v\ E X. 

The following lemma is a simple consequence of the minimality of the period 
chosen to represent an ultimately periodic set. 

Lemma 11. Let X C N be an ultimately periodic set of period px and preperiod 
ax- Let i,j > ax- If i ^ j mod px then there exists t < px such that either 
i + t S X and j + t ^ X or, i + t ^ X and j + t £ X. 

We assume that the reader is familiar with automata theory (see for instance 
[37] ) but let us recall some classical results. Let L C E* be a language over a finite 
alphabet S and i be a finite word over E. We set 

x~ l L = {z e E* | xz e L}. 

We can now define the Myhill-Nerode congruence. Let x,y S E*. We have x ~l V 
if and only if x~ 1 L = y~ 1 L. Moreover L is regular if and only if has a finite 
index being the number of states of the minimal automaton of L. 

Definition 12. For a sequence (Ui)i>o of integers, Nu(m) £ {1, . . . , to} denotes the 
number of values that are taken infinitely often by the sequence (Ui mod to)i>o- 

Proposition 13. Let U = (Ui)i>o be a positional numeration system satisfying 
©■ If X C N is an ultimately periodic U -recognizable set of period px, then any 
deterministic finite automaton accepting rep [/ (X) has at least Nu(px) states. 

Proof. Let ax be the preperiod of X. By Lemma [5l there exists L such that for 
any h > L, the words 

10>H ^Wl rep^W, i = 0,...,px-l 
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are greedy [/-representations. The sequence (Ui mod px)i>o takes infinitely often 
N := Nu(px) different values. Let hi, . . . , hx > L be such that 

i^ j ^ U hi ^ U h] mod px 

and hi, ... , hx can be chosen such that E7/, 4 > ax for all i G {1, ... , N}. 

By LemmalTT| for all i, j G {1, . . . , N} such that i ^ j, there exists t% j < px such 
that either U h% + t^j G X and fT/^ + $ X, or Uh t + Uj X and C/^. + £jj G X. 
Therefore, 

is a word such that either 

10i*-\T*Pu(t>x-V\ Witj g rep c/ (X) and io**-l ^^(pjic-i)!^. ^ rep£/ (X), 

or 

lO^-I^PcCf*- 1 )!^ £re Pc/ (X) and lO^'" 1 ^PtrCpx-i)!^ g rep v (X). 
Therefore the words lO^-\^Pu(px-i)\ )m _ _ ^ 10 ftjv-|re Pc7 (p x -i)| are pairwise non equi- 
valent for the relation r ^ TCPu (x) an d the minimal automaton of Tepjj(X) has at least 
N = Nu(p x ) states. □ 

The previous proposition has an immediate consequence for getting a bound on 
the period of a periodic set accepted by a given DFA. 

Corollary 14. Let U = (£/j)i>o be a positional numeration system satisfying (ffi). 
Assume that 

lim Nij(m) = +oo. 

m — >+oc 

Then the period of an ultimately periodic set ICfj such that Tepu(X) is accepted 
by a DFA with d states is bounded by the smallest integer sq such that for all m > sq, 
Nu(m) > d. 

A result similar to the previous corollary (in the sense that it permits to give an 
upper bound on the period) can be stated as follows. This result will be used later 
on in our decision procedure to compute some explicit estimate. 

Proposition 15. Let U = (Ui)i>o be a positional numeration system satisfying 
(PI) and X C N be an ultimately periodic U -recognizable set of period px- Let c 
be a divisor ofpx- If 1 occurs infinitely many times in (Ui mod c)i>o then any 
deterministic finite automaton accepting repu(X) has at least c states. 

Proof. Let ax be the preperiod of X. Applying several times Lemma [5j there exist 
ni, . . . ,n c such that 

10 n. 10 ne-i . . . 10 m l re P[r ( P x-l)|-|re Pc; (t)| repa ( t ) ) t = 0,...,p X -l 

are greedy [/-representations. Moreover, since 1 occurs infinitely many times in the 
sequence (Ui mod c)j>o, n\, . . . , n c can be chosen such that, for all j = 1, . . . , c, 

val[/(10™' ...10™ 1+ l re M px - 1 )l) = j mode 

and 

val C /(10" 1+|rcPt ' (px " 1)l ) > a x . 
For i,j g {1, . . . , c}, i ^ j, by Lemma [TT1 and since c divides px, the words 

10" 1 • • • 10™ 1 and • • • 10" 1 

are nonequivalent for ~ rep „(x)' This can be shown by concatenating some word 
of the kind fjl te Pubx-i)\-\ re Pt7 (t)| rePt7 (t) with t < p x , as in the proof of Proposi- 
tion This concludes the proof. □ 
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Definition 16. For a sequence (Ui)i>o of integers, if (£7$ mod m)j>o is ultimately 
periodic, we denote its (minimal) preperiod by (.[/(m) (we choose the notation i 
to allude to the word index which is equally used as preperiod) and its (minimal) 
period by 7rj/(m). 

Remark 17. Observe that for any linear recurrence sequence of order k satisfying 
(HJ), we have 

Nu{m) < 7Tu(m) < (Nu(m)) k . 
Therefore, lim m ^ +00 Njj(m) = +oo if and only if lim m _> +00 wu(rn) = +oo. Notice 
that if to = p.q with gcd(p, q) = 1, then 7T[/(to) = lcm{7T[/(p), iru(q)}. 

Now we want to obtain an upper bound on the preperiod of any ultimately 
periodic [/-recognizable set recognized by a given DFA. 

Proposition 18. Let U = (t/j)j>o be a linear numeration system. Let X C N 
be an ultimately periodic U -recognizable set of period px and preperiod ax- Then 
any deterministic finite automaton accepting repu(X) has at least \ repjj(ax — 1)| — 
i>u(j?x) states. 

The arguments of the following proof are similar to the one found in |26j . 

Proof. W.l.o.g. we can assume that \ repu(ax — 1)| — iu(px) > 0. The sequence 
(Ui mod px)i>o is ultimately periodic with preperiod ijjipx) and period iru(px)- 
Proceed by contradiction and assume that A is a deterministic finite automaton 
with less than | rep (7 (aj5f — 1)| — iuipx) states accepting repjj^X). There exist words 
u>, W4 such that the greedy [/-representation of ax — 1 can be factorized as 

rep;y(ax — 1) = ww^ 

with \w\ = |rep (7 (ax — 1)| — iu(px)- By the pumping lemma, w can be written 
W1W2W3 with W2 £ and for all i > 0, 

w\w 2 w^wa € rep (7 (X) W1W2W3W4, £ vep u (X). 

By minimality of ax and px , either ax -lei and for all n > 1 , ax + npx — lf-X, 
or ax — l&'X and for all n > 1, ax + npx — lEX. Using the ultimate periodicity 
of (Ui mod px)i>o, we observe that, because |u>4| = iuipx), for all i > 0, we have 

va\ u (w 1 wl 7ru[px) w 2 w 3 w 4 ) ee valuiw^wsw^+ivaluiw^^O^ 2 ™ 3 ^) mod p x . 

Therefore repeating a factor of length multiple of itu{Px) exactly px times does 
not change the value mod px and we get 

valu (wiw^* ™ u ^ Px ^ W2W3W4) ee valu (W1W2W3W4) mod px, 
leading to a contradiction. □ 

For the sake of completeness, we restate some well-known properties of ultimately 
periodic sets (see for instance [37] for a prologue on Pascal's machine for 6-ary 
systems) . 

Lemma 19. Let a, 6 be nonnegative integers and U = (C/j)j>o be a linear numera- 
tion system. The language 

valy 1 (aN + b) = {w £ A* v \ vaby(u;) GoN + dJcip 

is regular. In particular, ifN is U -recognizable then a DFA accepting repy(aN + 6) 
can be obtained efficiently and any ultimately periodic set is U -recognizable. 

Before giving the proof, notice that for any integer n > 0, val^ 1 (n) \ + A^j is a 
finite set of words {xi, . . . , Xt n } over Ajj such that val[/(xj) — n for alH = 1, . . . , t n . 
This set contains in particular rep[/(n). 
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Proof. Since regular sets are stable under finite modification (i.e., adding or re- 
moving a finite number of words in the language), we can assume that < b < a. 
The sequence (Ui mod a)i>o is ultimately periodic with prcpcriod £ = Lu(a) and 
period p = iru(a). It is an easy exercise to build a deterministic finite automaton A 
accepting the reversal of the words in {w S A^ \ v&\u(w) £ aN + b}. The alphabet 
of the automaton is A\j. States are pairs (r, s) where < r < a and < s < £ + p. 
The initial state is (0, 0). Final states are the ones with the first component equal 
to b. Transitions are defined as follows 

Vs < £ + p — 1 : (r, s) (jU s + r mod a, s + 1) 

(r,£ + p — l)-?-+(jUi + p-x + r mod a, £), 

for all j € Ajj. Notice that A does not check the greediness of the accepted words, 
the construction only relies on the [/-numerical value of the words modulo a. 

For the particular case, one has to consider the intersection of two regular lan- 
guages repy (N)nval^ 1 (aN + &). □ 

Remark 20. In the previous statement, the assumption about the [/-recognizability 
of N is of particular interest. Indeed, it is well-known that for an arbitrary linear nu- 
meration system, N is in general not [/-recognizable. If N is [/-recognizable, then U 
satisfies a linear recurrence relation [39], but the converse does not hold. Sufficient 
conditions on the recurrence relation that U satisfies for N to be [/-recognizable are 
given in [25j . 

Our decision procedure will also make use of the following result. 

Lemma 21. Let U = (C/i)i>o be an increasing sequence satisfying a linear recur- 
rence relation of order k of the kind ([TJ). The following assertions are equivalent: 

(i) liirim-j+oo Nu(m) = +oo 

(ii) for all prime divisors p of ak, lim,,->.-|-oo Njj(p v ) = +oo. 
In particular, if ak = ±1, then \im m ^ +oc Njj(m) = +oo. 

Proof. It is enough to show that (ii) implies (i). Let the prime decomposition of 
|ofc| be |dfc| = p^ 1 ■ ■ -Pr r with u\, . . . ,u r > 0. It is obvious that if m = p^ 1 ■ ■ -p v r r c 
with V\, . . . ,v r > and gcd(a^, c) = 1 then 

nu(m) = lcmjvrc/^ 1 ), . . . , 7Tt/(p" r ), 7T[/(c)}. 

Notice that m tends to infinity if and only if at least one of the Vj's or c tends to 
infinity. 

Assume first that for some j G {1, . . . , r}, Vj — > +oo. By assumption we have 
that lim^.^+oo Nu(p V j J ) = +oo and by Remark ll71 we get that lim^.^+oo ^uiPj 3 ) — 
+oo. Therefore 7rj/(m) takes values larger than any constant by considering an inte- 
ger m which is divisible by a sufficiently large power of pj. Again using Remark 1 171 
the same conclusion holds for Njj{m). 

Let C — {cq < c\ < C2 < ■ • • } be the set of natural numbers prime to ak- 
For all c S C, the sequence (Ui mod c)i>o is ultimately periodic but it is even 
purely periodic. Indeed, for all i > 0, Ui+k is determined by the k previous terms 
[/j+fc-i, ■••,[/«• But since gcd(afe, c) = 1, ak is invertible modulo c and for all i > 0, 
Ui mod c is also determined by the k following terms E/j+i, . . . , t/j+fc. By definition 
of Nu(c), the sequence (Ui mod c)i>o takes exactly Nu(c) different values because 
any term appears infinitely often. Let a be the function mapping m € N onto 
the smallest index a(m) such that U a r m -\ > m. Since U is increasing, a is non- 
decreasing and lim m ^ +00 a(m) — +oo. From this last observation and from the 
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pure periodicity of (Ui mod c)i>o, it follows that for all c € C, Nu(c) > a(c) 
because Uq < ■ ■ ■ < ?7 Q ( c )_i < c. Consequently, we get 

lim Nu(c n ) = +00. 

n — ' + 00 

Any large enough integer m contains either a large power of some pj or some 
large c prime to and consequently (i) holds. □ 

Theorem 22. Let U — (Ui)i>o be a linear numeration system such that N is 
U -recognizable, satisfying condition @. Assume that 

lim Nu(m) = +00. 

Then it is decidable whether or not a U -recognizable set is ultimately periodic. 

Proof. The sequence U satisfies a recurrence relation of order k of the kind |T]) . Let 
the prime decomposition of \ak\ be \a^\ = p" 1 • ■ -p" r with u±, . . . , u r > 0. Consider 
a DFA A with d states accepting a [/-recognizable set ICR Assume that X is 
ultimately periodic with a period 

PX =Pi ■■ -P r C 

where gcd(afc, c) = 1 and Vx, . . . , v r > 0. 

Since gcd(afc,c) = 1, with the same reasoning as in the proof of the previous 
lemma, the sequence (Ui mod c)i>o is purely periodic. Therefore, Uq = 1 appears 
infinitely often in (Ui mod c)i>o- Since c is a divisor of px, we can use Proposi- 
tion [15] and get c < d. 

By Proposition Q21 we get Njj(px) < d. Let j € {1, . . . , r}. Using Remark [T71 
we obtain 

Nu(p] 3 ) < 7Tu(p7) < *u(px) < (Nu( Px )) k < d k . 
The assumption lim„ w+00 Njj(m) = +00 implies that \\m v ^ +oa Nu(p V j) = +00. 
Observe that Njj(p'j) < Njj(p™) whenever v < w. Consequently the exponent Vj 
occurring in the decomposition of px is bounded by Sj where Sj is the smallest 
integer such that for all v > sj, Nu(Pj) > d k . This bound Sj can be effectively 
computed as follows. For any v, Njj(p^) can be found in a finite number of opera- 
tions by inspecting the first values of (Ui mod Pj)i>o and looking for two identical 
fc-tuples made of k consecutive elements. Once the period is determined, one imme- 
diately gets the values that are repeated infinitely often. Since the map v <— > Nu(p^) 
is non-decreasing, one has to compute Nu(pj) < Nu(p^) < ■ ■ ■ until finding the 
first value Sj such that Njj(pj 3 ) > d k . 

If X is ultimately periodic, then the admissible periods are bounded by the 
constant 

/• i'V ■■■/>: d 

which is effectively computable. Then, using Proposition [THl the admissible prepe- 
riods ax must satisfy 

\rep u (a x - 1)| < d + max^O)) 
p<p 

where |rep;j(a)| < | rep t/ (6) | whenever a < b. This last observation shows that a 
bound on the admissible preperiods of X can be given effectively. 

Consequently the sets of admissible preperiods and periods that we have to check 
are finite. For each pair (a,p) of admissible preperiods and periods, there are at 
most 2 a 2 p distinct ultimately periodic sets. Thanks to Lemma [T9l one can build 
an automaton for each of them and then compare the language L accepted by 
this automaton with repu(X). Recall that testing whether L \ rep [/ (A) = and 
repjj(X) \ L = is decidable algorithmically. □ 
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In view of the previous result, it is natural to characterize linear recurrence 
sequences U such that lim m ^ +00 Nu(m) = +00. It is exactly the aim of Section [31 
For instance, we have the immediate special case that was treated in [8]. 

Corollary 23. Let U = (L^)i>o be a linear numeration system such that N is 
U -recognizable and satisfying a recurrence relation of order k of the kind (|TJ) with 
ah = ±1 and condition It is decidable whether or not a U -recognizable set is 
ultimately periodic. 

Proof. It is enough to observe that if = ±1, then \im m ^ +oc Njj(m) = +00, see 
Lemma □ 

Example 24. Let a G (0, 1) be irrational, with simple continued fraction a = 
[1, <zi, a, 03, a, . . .], that is, 02% = a (i > 1), where a, an-\ € N \ {0} for all i > 1. 
The numerators of its even-indexed convergents are given by the recurrence XJii = 
(aa,2i-i + 2)L?2i-2 — C^2i-4 {i > 1) with initial conditions U-2 = 1 — a, Uq = 1. 
It was shown in |16j that every positive integer n has a unique representation of 
the form n = ^2 i>0 dufJn, where the digits satisfy < c?2i < cio-2i+i + 1, and the 
additional condition: If for some < k < £ the digits d2k and &2t attain their 
maximum values, then there exists j satisfying k < j < £ such that d%j < ao2j+i- 
It was recently employed in [3Hj, where it was dubbed ^-sequences. The special 
case a — <Z2j_i — 1 for all i > 1 is an exotic ternary numeration system since then 
< c?2i < 2. The additional condition then states that between any two digits 
2 there must be a digit 0. Also a = (1 + \/5)/2 is the golden ratio and the C/^ 
are the even- indexed Fibonacci numbers. This special case was used to investigate 
irregularities of distribution of sequences [TU], [TT]. In [TJ] it was used to compress 
sparse bit-strings, in [TT] for answering a combinatorial question raised in [B], and 
in [18] for providing a polynomial-time algorithm for combinatorial games. For all 
these numeration systems Problem[Tjis decidable, since the coefficient of L*2i-4 is —1. 
The same conclusion holds for the numeration system based on the denominators 
of the even convergents, also given in [16) . 

Remark 25. We have thus obtained a decision procedure for our Problem [T] when 
linim—^oo Njj (m) = +00 and in particular when the coefficient occurring in |T]) 
is equal to ±1. On the other hand, whenever gcd(ai, . . . , a^) = g > 2, for all n > 1 
and for all i large enough, we have Ui = mod g n and the assumption about 
Njj(m) in Theorem 1 2 2 1 does not hold. Indeed, the only value taken infinitely often 
by the sequence (Ui mod g n )i>o is 0, so Njj(m) equals 1 for infinitely many values 
of m. Notice in particular, that the same observation can be made for the usual 
&-ary system (b > 2) where the only value taken infinitely often by the sequence (b l 
mod b n )i>o is 0, for all n > 1. 

To conclude this section, we make a small digression. We show how to use a result 
of Engstrom about preperiods [14] to get some special linear numeration systems 
where lim m _> +00 Njj(m) = +00. In Engstrom's paper the problem of finding a 
general period for a given recurrence relation modulo m for any initial conditions 
is considered. Notice that in 4DJ, M. Ward considers the problem where the initial 
conditions are fixed and then the period modulo m has to be determined. 

Theorem 26. 14, Theorem 9] Let U = (Ui)i>o be a linear recurrence seguence of 
order k of the kind ([T]) and p be a prime divisor of a^. If there exists s(p) < k such 
that dfc, . . . , afc_ s ( p ) +1 = (mod p) and ak- s (p) ^ (mod p), then Lu(p v ) < vs(p). 

Remark 27. Assume that we are dealing with a linear numeration system U = 
(Ui)i>o satisfying fTJ and that the assumptions of the previous theorem hold for 
all prime divisors p of ctfc (which is equivalent to the fact that gcd(ai, . . . , a^) = 1). 
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Let xu{x) be the characteristic polynomial of U as defined in Q. Assume that 
f3 > 1 is a root of multiplicity £ > 1 of Y[/(x) satisfying: 

• for any other root 7 £ C of I7I < /3, 

• (3 < p 1 /^) for all prime divisors p of a^. 

There exists some constant c such that Ui ~ ci l_1 /3\ Let p be a prime divisor 
of afc and j p (w) be the largest index j such that Uj < p v . Let t > s(p) be a 
real number such that (3 < p 1 ! 1 < p 1 / 8 ^). For v large enough, we have U\ vt \ < p° . 
Consequently, for v large enough, j p (v) > J . From the previous theorem, we have 
i'u(p v ) < vs(p), Therefore for v large enough, U Lu ^) < ■■■ < Uj ^ are the first 
terms of the periodic part of (Ui mod p v )i>o and Njj(p v ) > L W ^J — ws(p) + 1. This 
means that for all prime divisors p of a^, Njj(p v ) — > +00 as v — ► +00. Therefore, 
by Lemma lim m ^ +00 Njj(m) = +00 and we can apply our decision procedure 
given by Theorem [22] (whenever N is {/-recognizable) . 

Example 28. Consider the linear recurrence sequence given by t/,+3 = Ui+i + 3Ui 
for i > and Ui = i + 1 for i = 0, 1, 2. The first terms of the sequence are 



With the above notation, (3 ~ 1.6717 < 3 and the other two complex roots have 
modulus close to 1.34. We also have s(3) = 1. Thanks to Theoremf^!! the preperiod 
i<u{3 v ) is bounded by v. On the other hand, we have Ui ~ c(3 l for some c > 0. Notice 
that /3 < 3 1 / 2 < 3. Therefore for v large enough, L^-y ~ c/3 2u < 3". Consequently, 
the elements U v < ■ ■ ■ < U% v appear in the periodic part. In the following table, 
these elements have been underlined. 



V 


preperiod 


period 


3 


1,2,3 


(5, 9, 14, 24, 14, 12, 5, 0, 14, 15, 14, 3, 5, 18, 14, 6, 14, 21) 


4 


1,2,3,5 


(9, 14, 24, 41, 66, 32, 27, 68, 42, 68, 3, 32, 45, 41, 60, 14, . . .) 


5 


1,2,3,5,9 


(14, 24, 41, 66, 113, 189, 68, 42, 149, 3, 32, 207, 41, 60, 176, . . .) 



3. Linear Recurrence Sequences and Residue Classes 

As was observed in Remarkl25[ since our approach to solve the decision problem 
is requiring that lim m ^ +oc Njj(m) = +00, it can only be applied to linear recur- 
rence sequences (TTJ) for which gcd(ai, . . . ,Ofc) = 1. In this section, our aim is to de- 
termine which linear recurrence sequences U are such that lim„ w+00 Njj(m) = +00. 
To that end, it is clear (see Lemma l2"Tj) that we have only to focus on the behavior 
of Njj(p v ) for any prime p dividing afc. 

Throughout this section we let (£/j)j>o be a linear recurrence sequence satisfying 
(fTJ). We assume that (Ui)i>o satisfies no recurrence of smaller order than k. It is 
well-known (this result is sometimes referred as Kronecker's theorem, see [5S]) that 
under the assumption afc ^ 0, this is equivalent to assume that 



Let p be a prime number. We recap some background on the p-adic numbers 
(see for instance [H]). We can put an absolute value | • | p on Z as follows. For 
each integer n ^ 0, we can write n = p v £ with I such that gcd(p, i) = 1. We define 
\n\ p = p~ v and |0| p = 0. We note that this absolute value extends on Q by declaring 
\a/b\ p = \a\p/\b\ p for a, b e Z, b ^ 0. In particular, for all a, b £ Q, |a.6| p = |o| p .|6| p . 
Note that this absolute value is non-Archimedian; that is, it satisfies for all a, b £ Q 



1,2,3, 5, 9, 14, 24, 41, 66, 113, 189, 311, 528, 878, 1461, 2462, 4095, 



det 









\a + b\ p < max{|a| p , |6| p }. 
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If we complete Q with respect to this absolute value, we obtain the field of p- 
adic rationals, which we denote by Q p . We can take the algebraic closure of Q p ; 
the absolute value | • \ p extends to this algebraic closure. The algebraic closure is 
not complete, however. Completing this algebraic closure, we obtain a complete 
algebraically closed field C p with absolute value | • | p , which restricts to the p-adic 
absolute value on Q. The closed unit ball 

{x € Q P : \x\ p < 1} 

is called the set of p-adic integers and we denote it by Z p . The ordinary integers 
are dense in Z p . The p-adic rationals can be viewed as the formal expressions of 
the form 

c-nP~ N H 1- c-ip^ 1 + c + cip + c 2 p 2 H 

where Cj € {0, . . . ,p — 1}, N G Z. The p-adic integers are identified with formal 
expressions involving only non-negative powers of p. 

We let 

(4) xu(x) = x k - axx k ~ l a k 

be the characteristic polynomial of (Ui)i>o and 

(5) Pu{x) = x k xu(l/x) = 1 - aix a k x k . 

For a k 7^ 0, observe that if a±, . . . , a s are the roots of xu, then the reciprocals 
1/cti, ■ • ■ , l/a s are exactly the roots of Pjj. 

Remark 29. If (£/j)j e t*j satisfies ([TJ) with initial conditions J7o, ■ ■ ■ , U k -i, then the 
ordinary power series generating function is a rational function: 

,fc-i 



The converse also holds. The sequence of Taylor coefficients of any rational function 
R(x)/Pjj(x) where R(x) is a polynomial of degree less than k, satisfies ([IJ. See for 
instance [T3 p. 6]. 

Our goal is to prove the following result. Notice that similar developments can 
be found in [34] . 

Theorem 30. We have Njj(p v ) y4 +oo as v — * +oo if and only if Pjj(x) — 
A(x)B(x) with A{x), B(x) € 7L\x\ such that: 

(i) B(x) = 1 (mod pZ[x\); 

(ii) A(x) has no repeated roots and all its roots are roots of unity. 
In that case, we have moreover A(0) = B(0) = 1. 

Proof. We note that one direction is fairly simple. Assume that Pjj{x) has such a 
factorization Pjj(x) — A{x)B{x). From (ii), there is a natural number d such that 
A(x) | (x d — 1). In view of Remark 1^1 there exist some polynomials Q(x), R(x) such 
that 

( X^TT * - (^-VQW _ (x d -l)Q(x) _ R(x)_ 

[ ^ Mx) A(x)B{x) B(xY 

By assumption (i), there exists some integer polynomial B\(x) such that B(x) = 
1 — pBi(x). Hence 



Ei=0 Ui X 1 J2i+j<k ai ^3 X 



(x«-l)^U iX * = ; = ^R(x) Bl (x) 

i>0 y y 1 i>0 
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In particular, for any fixed integer v, (x d ~ 1) X)i>o ^ x% ^ s congruent to a polynomial 
(mod p v ). This means that Ui+d = U (mod p") for all i sufficiently large. In 
particular, there are at most d values which can occur infinitely often mod p v ; that 
is, Nu(p v ) < d for every v. 

To do the other direction is a little more work and we use p-adic methods. We 
first note that v <— > Njj(p v ) is a non-decreasing function, i.e., 

(6) Nu(p w ) > Nu{p") whenever w > v. 

In particular, if N(p v ) -f± +oo then there is some d such that Nu(p v ) = d for all v 
sufficiently large. We can pick integers ai tV , . . . , a^ v such that if Ui = a (mod p v ) 
for infinitely many i, then a = a,j tV (mod p v ) for some j £ {1, . . . , d}. Since 

{(ai )TO modp v ), . . . , (a d , w modp 11 )} = {(ai,„ modp"), . . . , (a d , v modp 1 ')} 

for w > v, there is no loss of generality to assume that 

= dj> (mod p°) 

for w > v, 1 < j ; < d. It follows that for 1 < j < d the sequence (a,j }V ) v >i is Cauchy 
in Tip. Thus there exist b\, . . . , bd € Z p such that for 1 < j ; < rf, 

(7) a^t, — > 6j as in +oo. 
Let 

^ - n(ui-6,-)GZp. 

i'=l 

Note that since the set of linear recurrence sequences is closed under (Hadamard) 
product, sum, and scalar multiplication (see for instance |15j). the sequence (Vj),->o 
satisfies a linear recurrence over Z p . By construction, Vj is eventually in p v Z p for 
any fixed v, since any values of U that are not congruent to one of b\, . . . , 6<j mod 
can only occur finitely many times. That is for any v, \Vi\ p < for i sufficiently 
large. Hence |V^| P — > as « — > +00. Since (Vi)j>o satisfies a linear recurrence, the 
power series 

V(x) := Vrf 

i>0 

is a rational power series in Q p (x). Moreover, V(x) converges on the closed unit 
disc Z p , since \Vi\ p — > (in non-Archimedian fields this is enough to guarantee 
convergence: a series Yli>o~fi converges in Q p if and only if lim^ +OQ |7j| p = 0). 
Since V(x) is a rational series and it converges on the unit disc, its poles j3i, . . . , fi r € 
C p must satisfy \/3j\ p > 1, for I < j < r. 

To finish the proof, we will make use of Lemma I5T1 stated below. For the sake of 
clarity, we have separated this technical result from this proof. We note that the 
statement of this lemma is very close to what we already have, but it makes the 
additional assumption that the poles of U(x) generate a free Abelian subgroup. In 
general, the poles of U(a;) generate a finitely generated Abelian subgroup of C p . 
From the so-called fundamental theorem of finitely generated Abelian groups (see 
for instance p. 141]), this group is isomorphic to Z e x T, for some finite Abelian 
group T and integer e > 0. 

Let us show how to get rid of the torsion group T to be able to invoke Lemma CTT1 
Let a = #T. For < b < a, instead of taking the sequence (J7i)i>o, consider the 
sequence (£7^ )<>o := (U a i+b)i>a- This latter sequence satisfies a linear recurrence 
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and the poles of the generating function \J^(x) of (U^)i>o are the a'th powera^of 
the poles of U(x). Consequently, the poles of \j( b \x) generate a finitely generated 
torsion-free Abelian group, which is necessarily a free Abclian group. 

Since the poles /3i, . . . ,j3 r € C p of V(x) satisfy \/3j\ p > 1, for 1 < j < r, with the 
same reasoning, we get that the poles of the rational function 

i>0 

are a'th powers of . . . ,(3 r and for 1 < j < r, \(3j\ P = \Pj\ p > 1- We can invoke 

Lemma I3T1 applied to the sequence (U i ; )i>o (it satisfies some linear recurrence 
sequence with integer coefficients), and we deduce that any pole 7 of U^ b \x) satisfies 
either \j\ p > 1 or 7 = 1. In particular, the distinct poles 1/ai, . . . , l/a s of U(x), 
with a>i,...,a s S C p , being a'th roots of the poles of U^(x), are either such that 
|l/aj|p > 1 or roots of unity. By minimality of the order k of the recurrence 
satisfied by (Ui)i>a, the poles of U(a;) are precisely the roots of Pu{%)- We factor 

P v (x) = (1 - Six) ■ ■ • (1 - S k x) 

(each Sj is one of oti, . . . , a s , although they may be repeated). Let us factor Pjj{x) 
as A(x)B(x) where 

A(x) = Yl i 1 - 5 J X ) and B ( x ) = n t 1 ~ 
b':|*j|p=i} {r-\Sj\ P <i} 

By assumption Pjj{x) £ Z[x]. Moreover if if is a splitting field of Pu(x) over Q 
then any automorphism of K must permute the set of Sj with \Sj\ p < 1, since 
the automorphism permutes the entire set of Sj's and it must send roots of unity 
to roots of unity. Thus B{x) is a rational polynomial, since it is fixed by every 
automorphism of K. Note that if n > 0, then the coefficient of x n in B(x) is given 
by a sum of products of n elements in {Sj : \5j\ p < 1}. The set of algebraic integers 
is a subring of C p and the only rationals that are algebraic integers are in fact 
integers. Since the o^'s are algebraic integers, P>{x) is thus an integer polynomial. 
Moreover, since the p-adic absolute value is non-Archimedian, the coefficient of x n 
in B(x), n > 0, has p-adic absolute value strictly less than 1. Note that an integer 
m satisfying \m\ p < 1 is necessarily a multiple of p. Hence B(x) = 1 mod pZ[x]. 

Let us turn to the polynomial A(x). The roots of A(x) are roots of unity. 
Moreover, A(x) € Z[x] by the same reasoning as before. 

To finish the proof, we have one last thing to show: we need to know that A(x) 
has no repeated roots. To do this, we show that the poles of U(a;) that are roots of 
unity are simple. Recall that l/cei, . . . , l/a s are the distinct poles of U(a;). We may 
assume that there exists t > such that a\, . . . , at are roots of unity and that for 
j = t + 1, . . ., s, \otj\p < 1. Then there exist polynomials qj € C p [x], j = 1, . . . , s, 
such that we have for all i, 

s t s 

u i = q i w a ) = 12 q 3 (*) a ) + 12 * w a ) ■ 

i=i j=i f j=t+i ^ 

:=T< :=Wi 

Since for j = t + 1, . . . , s, \otj\ p < 1, we get that \Ui — Tj| p = \Wi\ p — > as i — > +00. 
Indeed, for any j, : i G N} is bounded by a constant. Since for 1 < j < t, 

^li we consider the exponential sum Ui = where the qj's are polynomials in 

Cj)[x] and the a^'s are the roots of the characteristic polynomial of (J7i)i>o, i- e -i the reciprocals 
of the poles of U(ir), then we obtain that = U ai+b = J2j=i Q^-PjO* + b) {oPf. Hence the 

poles of (x) are the a'th powers of the poles of U(x). 
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a.j is a root of unity, there exists a natural number a such that a" = 1 for all 
j € {1, . . . , *}. As before, we let T t (b) = T ai+b for < b < a. Thus 

t t 
T? ] = ^2q j (ai + b)af +b = £ a* (ai + 6) 

is a polynomial with coefhcients in C p denoted by 

Let e > 0. By definition of the frj's (see © above) for all large enough i, there 
exists € {1, . . . ,d} such that \Ui — iW)|p < e. Since \Ui — Tj| p — * as i — > +oo, 
we get 

|(Ti -&!)••• (T< - 6 d )|p < JJd^ - Cilp + \Ui - h{i)\v + \h{i] ~ bj\ P ) 

3=1 

as i — > +oo because every factor is bounded by a constant and one tends to zero. 
Thus for < b < a, as = gb(i), we have 

K5&CO -&].)••• (sfb(») - Ml? ° as * +°°- 

Consider the polynomial defined by 

h(i) ■= (56 (*) -&].)••• (fif5(i) - M- 

Let no be a natural number. Notice that by the Binomial Theorem, we get \hb(no + 
p") — hi(no)\p < Cp~ v — > as v — ► +00, where C is the maximum of the p-adic 
absolute values of the coefficients of hb- So since |/i&(z)| p — + as i — > +00, we get0 
that hb{no) = 0. Thus for all integer i each /i&(i) = 0, < b < a, and so each 
gb is a constant polynomial (this is a consequence of the Fundamental Theorem of 
Algebra, hb{i) = for all i implies that gb takes infinitely often the same value 
amongst b\, . . . ,bd so it is constant). It follows that Tj is periodic: 

\fi > 0, T l+a = Tj. 

So (x a — l)U(x) has no poles on the closed unit disc so in particular on the unit 
circle, as 

\U i+a - Ui\ p < \T t+a - Ti\ p + \W i+a - Wi\ p -> 
as i — » +00. We know that U(x) = Q(x)/Pu(x) with Q and P[/ relatively prime (by 
minimality assumption on order of recurrence satisfied by (£/j)j>o). Consequently, 
A(x) divides (x a — 1) because (x a — 1)Q(x)/Pjj(x) has no poles on the unit circle. 
This shows that it has no repeated roots, completing the proof. □ 

Let us consider the technical lemma used in the previous proof. 

Lemma 31. Let (Ui)i>o be an integer linear recurrence sequence satisfying 
and let U(x) = X^>o ^i x% ^ ^p[[ x ]] ^ e ^ e corresponding rational power series. 
Suppose that the multiplicative subgroup of C* generated by the (finitely many) 
poles of U(x) is a free Abelian group and let Vi — Yij=x(Ui ~ bj) for i > where 
b\, . . . , bet € lip. If the rational power series \/(x) = J2i>o ^ ^pIM] has poles 
f3%,...,f3 r € C p satisfying \/3j\ p > 1 for 1 < j < r, then every pole 7 S C p of U(x) 
either satisfies \^/\ p > 1 or 7 = 1. 

Proof. Let 1/ai, . . . , l/ct s be the distinct poles of U(x) with a\, . . . ,a s £ C p . Note 
that cannot be a pole of U. We first claim that \ctj\ p < 1 for 1 < j < s. To see 
this, notice that to be a pole, each 1/a.j must satisfy Pu(l/aj) = 0. This means 
that 

1 - axjoti - ■ ■ ■ - ak/cti = 0. 

2 Indced, \hf,(no)\ p < \hij(no) — h^no + p v )\ P + \hf,(no +p")|p and both terms tend to as 
V —* +00. 
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Consequently, for all j G {1, . . . , s}, we have 



\ai/aj-\ h a k /a k A 



1. 



By the non-Archimedian property we have \ai/a^\ p > 1 for some I G {1, . . . , k}. 

Since ag G Z, \ai\ p < 1. This gives the claim and so |oj| p < 1 



Thus \a 3 \ p < 



3 IP 

1 and |a,-|„ < 1 for 



for 1 < j < s. We may assume that |ai| p = • • • = \a t \ P — ± a,±±<a \^j\ p 
< t < j < s. So, we know that the poles of U have p-adic absolute value > 1. 
Next there exist polynomials qi(x), . . . , q s (x) G C p [x] such that 

s 

3 = 1 



Moreover, define 

Vi := (U, -&!)••• (U, - b d ) = c d Uf + c^Uf- 1 + • • • + c 
for some c , . . . , Cd-i, Cd = 1 G Z p . Hence by Multinomial Theorem, we have 



j=0 jiH \-3b=3 

3i,-,3s>0 



E 



- E 



3=0 



3l~\ \-3s=3 

3i,-,3s>0 




Since the roots of the characteristic polynomial are the reciprocals of the poles of 
the corresponding rational power series, it follows that the set of poles of V(x) is 
contained in the set 




ji, ■ ■ -,3s > 0, h 



+ js <d 



By assumption, the poles of V(x) all have p-adic absolute value strictly greater 



than 1. Note that 



IE 



-31 



-31 



n 



e=t+i ' 



-31 



is > 1 if and only 



if je > for some £ > t. Therefore, we can conclude that the possible poles of V 
supported only by products on a\, . . . , at do not occur. So for instance, af is not 
a pole of V. 

Let G denote the multiplicative subgroup of generated by a±, . . . , at- By 
assumption, G is a subgroup of a finitely generated free Abelian group, and hence 
G = Z e , for some natural number e > 0. To conclude the proof, it is enough to 
show that e = 0. Because in that case, one can conclude that the only possible pole 
7 of U such that \y\ p = 1 is 1. To see this, suppose that e > and let 71, . . . ,j e be 
generators for a free Abelian group of rank e. Then for 1 < i < t, we can write 



ai 



3=1 



where bij are integers. We relabel if necessary so that 

• |6 M | = max{|6 i4 |} > 0; 

• \bi t j\ = max{|6 i; j| : bu = b\ t t for I < j}. 

By construction, af cannot be written as a different word in a\, . . . , at of length 
at most d. Then the expression for Vi above has an occurrence of 

c dqi (i) d af 
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that cannot be canceled by any other pole, by our selection of ot\. Consequently, af 
should be a pole of V which contradicts the conclusion of the previous paragraph. 
So, e — and the result follows. □ 

As it was stressed in the introduction of this section, we are now able to check, 
whether or not the number of values taken infinitely often by a linear recurrence 
sequence modulo p v tends to infinity as v — > +00. Let us consider the following 
example. 

Example 32. Consider the linear recurrence sequence (Ui)i>o given in [5D] and 
defined by 1^+4 = 3t/,+3 + 2Ui + 2 + 3Ui for i > and C/j = i + 1 for i = 0, . . . , 3. As 
shown in |20j . addition within this linear numeration system is not computable by 
a finite automaton. Nevertheless, we can show by applying the previous theorem 
that Njj(3 v ) — > +00 as v — > +00. One has 

Pu{x) = 1 - 3a; - 2x 2 - 3x 4 

and it is not difficult to see that this polynomial cannot be factorized as A(x)B{x) 
with two factors satisfying hypotheses of Theorem [3U1 This example shows that 
our decision procedure given by Theorem [521 can take care of numeration systems 
not handled by 01301 El- 
Example 33. Consider the recurrence relation 

U t+5 = W l+i + 3U t+3 - U i+2 + 6U i+1 + Wi, Vi > 0. 
With the above notation, xu( x ) — x 5 — 6x 4 — 3x 3 + x 2 — Qx — 3 and 



Pu(x) = l-6x-3x 2 



6x 4 



3x 5 = (x 3 + 1) (-3x 2 - 6x + 1) . 



With the initial conditions Ui = i + 1 for i = 0, . 
does not satisfy any relation of shorter length as 





(I 


2 


3 


4 


5 \ 




2 


3 


4 


5 


54 


det 


3 


4 


5 


54 


359 




4 


5 


54 


359 


2344 




V5 


54 


359 


2344 


15129/ 



A(x) B(x) 

, 4, the corresponding sequence 



8458240 ^ 0. 



Even if the gcd of the coefficients of the recurrence is 1, 
assumptions of Theorem [30] for p = 3, Njj(3 v ) -/-* +00 as v 
table gives the first values of Nu(3 v ). 



since P\j satisfies the 
— > +00. The following 



V 


period 


^(3°) 


1 


(1,0,1,2,0,2) 


3 


2 


(4,0,1,5,0,8) 


5 


3 


(22,9,19,5,18,8) 


6 


4 


(49,63,19,32,18,62) 


6 


5 


(211,225,19,32,18,224) 


6 



4. A Decision Procedure for a Class of Abstract Numeration 

Systems 

An abstract numeration system S = (L, S, <) is given by an infinite regular 
language L over a totally ordered alphabet (S, <) [29.. By enumerating the words 
of L in genealogical order, we get a one-to-one correspondence denoted rep s between 
N and L. In particular, is represented by the first word in L. The reciprocal map 
associating a word w G L to its index in the genealogically ordered language L is 
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denoted vals. A set X C N of integers is S -recognizable if the language rep s (X) 
over £ is regular. 

Let S = (L, £, <) be an abstract numeration system built over an infinite reg- 
ular language L having Ml = (Qi> <Zo,£> S, #l, F£) as minimal automaton. The 
transition function 5l '■ Ql x E — > Qz, is extended on Q/, x S* by 5l(q,s) = q 
and <5i(g,aw) = <5l(<5l(<7, a), w) for all g £ Ql, a € S and w <E S*. We denote 
by Ui(q) (resp. Vj(q)) the number of words of length i (resp. < i) accepted from 
q £ Ql in Ml- By classical arguments, the sequences (uj(g))j>o (resp. (vi(q))i>o) 
satisfy the same homogenous linear recurrence relation for all q £ Ql (for details, 
see Remark [HS]). 

In this section, we consider, with some extra hypothesis on the abstract numer- 
ation system, the following decidability question analogous to Problem [T] 

Problem 2. Given an abstract numeration system S and a set X C N such that 
rep s (X) is recognized by a (deterministic) finite automaton, is it decidable whether 
or not X is ultimately periodic, i.e., whether or not X is a finite union of arithmetic 
progressions ? 

Abstract numeration systems are a generalization of positional numeration sys- 
tems U = (Ui)i>o for which N is [/-recognizable. 

Example 34. Take the language L — {e} U 1{0, 01}* and assume < 1. Ordering 
the words of L in genealogical order: e, 1, 10, 100, 101, 1000, 1001, . . . gives back the 
Fibonacci system. 

The following example shows that the class of abstract numeration systems is 
strictly bigger than the class of linear numeration systems for which N is recogniz- 
able. 

Example 35. Consider the language L = {e} U {a, ab}* U {c, cd}* and the ordering 
a < b < c < d of the alphabet. If we order the first words in L we get 






e 


5 


cc 


10 


ccc 


15 


aaba 


20 


cede 


1 


a 


6 


cd 


11 


ccd 


16 


abaa 


21 


edec 


2 


c 


7 


aaa 


12 


cdc 


17 


abab 


22 


eded 


3 


aa 


8 


aab 


13 


aaaa 


18 


cccc 


23 


aaaaa 


4 


ab 


9 


aba 


14 


aaab 


19 


cccd 


24 


aaaab 



Notice that there is no bijection v : {a, b, c,d} — > N between {a, b, c, d} and a 
set of integers leading to a positional linear numeration system. Otherwise stated, 
a,b,c,d cannot be identified with usual "digits". Indeed, assume that there exists 
a sequence U of integers such that for all x\ ■ ■ ■ x n G L, with Xi in {a, 6, c, d} for 
all i, v&\u(i/(xi) ■ ■ ■ v(x n )) = vsls(x\ ■ ■ ■ x n ). Since valg(a) = 1 and vals(c) = 2 we 
get Uq = 1, v{a) = 1 and v{c) = 2. Moreover, vals(aa) = 3 = v(a)Ui + u(a)Uo, 
so Ui = 2. Therefore, valy (^(c)^(c)) = 2U\ + 2Uq = 6 but vals(cc) = 5 getting a 
contradiction. 




Figure 1. A DFA accepting L. 
For alii > 1, we have u^qo.l) = 2Fi and u (<7o,l) = 1. Consequently, for i > 1, 

i i 

v t (qo,L) = l + ^2u n (q Q . L ) = 1 + 2^F„. 

n— 1 n—1 
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Notice that for i > 1, Vi(q 0:L ) - Vj_i(g ,z) = Uifao,!,) = 2F t . Consequently, by 
definition of the Fibonacci sequence, we get for all i > 3, 

Vi(<?o,i) ~ v<_i(go,i) = (v 4 -i((? ,l) - v;_ 2 (go,L)) + (v<_2(go,i) - v;_ 3 (q ,L)) 
and 

Vi(go,i) = 2vi_i(9o,L) - Vi_ 3 (g 0: L), with v (g ,L) = l,Vi(g ,i) = 3, v 2 (cfo,i,) = 7. 

Remark 36. The computation given in the previous example to obtain a homoge- 
nous linear recurrence relation for the sequence (vj(go,i))i>o can be carried on in 
general. Let q G Ql- The sequence (ui(?))i>o satisfies a homogenous linear re- 
currence relation of order k whose characteristic polynomial is the characteristic 
polynomial of the adjacency matrix of Ml- There exist eti, . . . € Z such that 
for all i > 0, 

Uj+fe(?) = aiUj+fe_i(g) H ha fc Uj(g). 

Consequently, we have for all i > 

Vi +fe+ i(q)-Vi +fe (g) = u i+fe+ i(q) = ai(v i+ k(q)-v i+ k-i(q))-\ hafc(v i+ i(g)-Vj (<?)). 

Therefore the sequence (v,(q))j>o satisfies a homogenous linear recurrence relation 
of order k + 1 . 

As shown by the following lemma, in an abstract numeration system, the different 
sequences (Uj(g))j>o, for q £ Ql, are replacing the single sequence (C/j)j>o defining 
a positional numeration system as in Definition [T] 

Lemma 37. [25] £ef w = <7i • ■ ■ o n € L. We have 

\w\ 

(8) val<j(iu) = ^2PqA w ) u H~t(q) 
w/iere 

(9) Ai.iM : = #{°" < CT i I 8 L (qo,L, <J\ ■ ■ ■ Vi-iv) = q} + l g ,? ,i 
/or i = 1, . . . , |io| . 

Recall that lq, q > is equal to 1 if g = q' and it is equal to otherwise. 

Proposition 38. |29j Let S = [L, E, <) be an abstract numeration system built 
over an infinite regular language L over E. Any ultimately periodic set X is S- 
recognizable and a DFA accepting iej> s (X) can be effectively obtained. 

Recall that an automaton is trim if it is accessible and coaccessible (each state 
can be reached from the initial state and from each state, one can reach a final 
state). 

Proposition 39. Let S = {L, E, <) be an abstract numeration system such that 
for all states q of the trim minimal automaton Ml = (Ql, qo,L, E, 5l, Fl) of L, 

lim Ui(q) = +oo 

i — >+oc 

and Uj(go,i) > for all i > 0. If X C N is an ultimately periodic set of pe- 
riod px, then any deterministic finite automaton accepting rep s (X) has at least 
\N v (px)/#Ql~] states where v = (v;(g ,i,))i>o- 

Proof. Let ax be the preperiod of X. Since for all states q of Ml, we have 
lirm^ +00 Ui(q) = +oo, there exists a minimal constant J > such that Uj(g) > px 
for all q G Ql- Consider for any i > 0, the word 

Wi = rep s (vi(g ,L)), 



20 



J. BELL, E. CHARLIER, A. S. FRAENKEL, AND M. RIGO 



corresponding to the first word of length i + 1 in the genealogically ordered language 
L. Consequently, for i > J — 1, Wi is factorized as Wi — aibi with = J and we 
define qi := <5l(<7o,l, a i)- Notice that bi is the smallest word of length J accepted 
from qi. By definition of J, from each qi, there are at least px words of length J 
leading to a final state. If we order them by genealogical ordering, we denote the 
px first of them by 

bi = b ifi < b it i < ■■■ < bi iPx -i. 
Notice that for t £ {0, . . . ,px ~ 1}, we have 

vals(a l 6 ij t) = val s (a i 6 l ) + t = Vj(g ,i) + t- 

The sequence (vi(go.i) mod px)i>o is ultimately periodic and takes infinitely 
often N v (px) =■ N different values. Let hi, . . . , /ijy > J — 1 such that 

i 7^ 3 => («?o,l) ^ (q 0jL ) mod p x 

and for all i G {1, . . . , N}, V/j 4 (qo,z) > aj. We have 

rep s (v fcj (go,z,)) = Wh; = a^&fcj and q hi = S L (q 0:L ,a hi ). 

The elements in the set {qh t , ■ ■ ■ , Qh N } can take only #Ql different values. So at 
least a := \N/#Ql\ of them are the same. For the sake of simplicity, assume that 
they are qh 1 , ■ ■ ■ , ■ Consequently, for i, j € {1, . . . , &} and for all t = 0, . . . ,px~ lj 
we have = 6/^,*. For all i, j € {1, . . . , a} such that i ^ j, by Lemma [TT| there 
exists t i: j < px such that either v h .(q _ L ) + tij & X and Vh s (Qo,L) + tjj ^ -X" or, 
v hi(9o,L) + ^ -X" and v hj (q 0tL ) + t,j G X. Therefore, the words and a hj 
do not belong to the same equivalence class for the relation ~ re p s (x) ■ This can be 
shown by concatenating the word bh^ti = .u ■• Hence the minimal automaton 
of rep s (A) has at least a states. □ 

Corollary 40. Let S = (L, S, <) be an abstract numeration system having the 
same properties as in Proposition \39l Assume that the sequence v = (vj(go,i))j>o 
is such that 

lim N v (m) = +oo. 

m — > + oo 

Then the period of an ultimately periodic set ICN such that rep s (A) is accepted 
by a DFA with d states is bounded by the smallest integer sq such that for all m > so, 
N v (m) > d^fQr,, where Ql is the set of states of the (trim) minimal automaton 
ofL. 

Proposition 41. Let S = (L, S, <) be an abstract numeration system. If X CN 
is an ultimately periodic set of period px such that rep s (X) is accepted by a DFA 
with d states, then the preperiod ax of X is bounded by an effectively computable 
constant C depending only on d and px ■ 

Proof. Let A = (Q,qo,Y^,S, F) be a DFA with d states accepting rep s (X). As 
usual, J\4l — (Ql, Qo,l, S, 5l, F£) is the minimal automaton of L and for any state 
q € Ql, u,(g) is the number of words of length i accepted from q in A4l- Since 
(u,(g))j>o satisfies a linear recurrence relation, the sequences (u,-(g) mod px)»>o 
are ultimately periodic for all q £ Ql- As usual, we denote by i u (q)(px) (resp. 
7r u(q)(px)) the preperiod (resp. the period) of (uj(g) mod px)i>o- We set 

I(p x ) := max L u{q) (p x ) 



and 



P(Px) ■= lcm gSQi 7T u ( 9 )(px)- 
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For ax large enough, we have |rep s (ax — 1)| > d#Q L . By the pumping lemma 
applied to the product automatory Ax. M.l, there exist x, y, z with y ^ e, \xy\ < 
d#Q L , 5(q ,x) = 5(q ,xy), S L (qa tL ,x) = S L (q 0:L ,xy) and such that 

rep s (a x - 1) = xyz 

and for all n > 0, 

(10) xy n z e rep 5 (A). 

Since \xy\ is bounded by a constant, we also have \z\ > I(px) if a>x is chosen large 
enough. 

Since \z\ > I(px), using (J3]), ([9|) and for all q G Ql the periodicity of the 
sequences (u,(g) mod px)i>o, we have for all £ > that 

(11) vals(xy ipxP ^ px 'yz) = va\s(xyz) mod px. 

Let us give some extra details on how we derive identity (fTT|) . Assume x = x\ ■ ■ ■ x r , 
y = yi ■ ■ ■ y s and z = z\ ■ ■ ■ z t . For all n > 1, using (jSJl for w = xy n z, we get 
\w\ = r + ns + t and 



va\ s (xy n z) = \ ^2f3q,i( w ) u \w\-i(q) 



q&Q L \*=1 

+ X] Pq,i( w ) u \w\-i(<l) H 1" X! / ? g^( w ) u |«)|- 4 ('?) 

i=r+l i=r+(n— l)s+l 

r+ns+t \ 

i=r+ns+l / 

where the first (resp. second, third) line corresponds, as explained below, to the 
contribution of x (resp. y n , z). By definition ^ of the coefficients (3 qi i(w), we know 
that P Q: i(w) depends only on x\ but (3 q .2(w) depends only on X2 and on 6l(%,Li xi). 
Continuing this way, (3 q ^{w) depends only on x r and on 5 L (qo.L, %i • • • 2V- 1) an d f° r 
1 < j < s, P qtr+ j(w) depends on yj and on 8 L (q 0iL , zyi ■ • • Vj-i)- Now /3q )r+s+ i(u>) 
depends only on j/i and on S L (q ,L,xyi ■■■y s ) = S L (q QtL ,xy) = 5 L (q 0t L,x). This 
implies that P q , r+s+ j(w) — (3 q ^ r+ j{w) for all q € Ql and all j € {1, . . . , s}. This 
argument can be repeated with every copy of y appearing in w. Consequently, the 
previous expansion becomes 

/ r r+s n— 1 

val s (xy n z) = ^ ( ^2 &qA w ) u H-i(q) + ^ PqA w ) ^2 u H-i-js(q) 

Q&Ql K »=1 i=r+l j=0 

(*) 



X] / 3 9,i( w ) U |«.|- l (<?))- 
-r+ns+1 



Assume now that n = l+£pxP(px), with I > 0. For g S Ql and i = r+1, . . . , r + s, 
we have 

n-l lpxP(px) 

(*) = X! u M-»-j«(90 = + X! u h|-i-js(9) 

i=0 j=l 



^The automaton A X A4f, is defined as follows. For any state (g, q') in the set of states 
Q X Ql, when reading a 6 S, one reaches in .A X A4l the state (S(q, a), (5^ (q', a)). The initial 
state is (qo; ?o,l) and the set of final states is F X F^. Roughly speaking, the product automaton 
mimics the behavior of both automata A and M^. 
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and the second term is congruent to modulo px due to the periodicity of the 
sequences (ui(q) mod px)i>o (recall that in the case we are considering, \z\ = t > 
I{px))- Consequently, for n = 1 + £pxP(px), we have 

(r r+s 
^fcH u i«'i-«(«)+ X! 0qA w ) n \w\-i(q) 
»=1 i=r+l 
r+ns+i v 

+ X! P<i,i( w ) u \w\-i(<l)) modpx- 

It is then easy to derive (JUJ. 

We now use the minimality of ax to get a contradiction. Assume that ax — 1 is 
in X (the case not in X is similar). Therefore for all n > 1, ax + npx — 1 is not in 
X. From Q3JJ), for f >0we get xy lpxPi - Px) yz g rep s (X), but from |Qj} this word 
represents a number of the kind ax + fipx — 1 with n > which cannot belong to 
X. 

Notice that C can be effectively estimated as follows. One has to choose a 
constant C such that ax > C implies | rep s (ax — 1)| — c?#Ql > I(px)- Since the 
abstract numeration system S, the period px and the number d of states are given, 
I(px) and repg(n) for all n > can be effectively computed. □ 

Theorem 42. Let S = (L,T*, <) be an abstract numeration system such that for 
all states q of the trim minimal automaton Ml = (Ql, Qo,l, S, Sl, Fl) of L 

lim Ui(q) = +oo 

i — >oo 

and Ui(go,i) > for all i > 0. Assume moreover that v = (v^ (qo t L))i>o * s such 
that lim m ^ +oa N v (m) = +oo. It is decidable whether or not an S -recognizable set 
is ultimately periodic. 

Proof. The proof is essentially the same as the one of Theorem First notice 
that the sequence v = (vi(qo.L))i>o satisfies a linear recurrence relation of the kind 
(fTJ) having a^ as last coefficient. Moreover it is increasing, since Uj(go,L) > for all 
i > 0. 

Let the prime decomposition of |dfc| be |dfc| = p" 1 • • -p^ r with u%, ... ,u r > 0. 
Consider a DFA A with d states accepting an S- recognizable set Icff. Assume 
that X is periodic with period 

PX =Pi ---Pr C 

where gcd(afc, c) = 1 and vx, ■ ■ ■ , v r > 0. 

By Proposition [3^1 we get N v (px) < d#<5i- Using Remark [T71 we obtain 

N v {c) < ttv(c) < n v (p x ) < (N v ( Px )) k < (d#Q L ) k . 

Let a(m) be defined as the largest index i such that Vi(qo,L) < m. Notice that 
since the sequence v is increasing, the map m i— ► a(m) is non-decreasing and 
lim m ^ +00 a(m) = +oo. Since gcd(a/c,c) = 1, the sequence (Vj(go,i) mod c)i>o is 
purely periodic and N v (c) > a(c). Therefore, a(c) < (<i#<5i,)' £ and we can give 
effectively an upper bound on c. 

Now we can give upper bound on the Vj's. The assumption lim m ^ +oc N v (m) = 
+oo implies that lim^^+oo N v (p^) — +oo and we have exactly the same reasoning 
as in the proof of Theorem l22l 

We have shown that if X is ultimately periodic, then its period px is bounded 
by a constant that can be effectively estimated. Using PropositiondTJ its preperiod 
is bounded by a constant which can also be computed effectively. 
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Consequently, the sets of admissible periods and preperiods we have to check are 
finite. Thanks to Proposition [33 one has to build an automaton for each ultimately 
periodic set corresponding to a pair of admissible preperiods and periods and then 
compare the accepted language with rep s (X). □ 

Example 43. The abstract numeration system given in Example [35] satisfies all 
the assumptions of the previous theorem. 

Theorem 0U can be used to decide particular instances of the HDOL periodicity 
problem. Let A, T be two finite alphabets. Consider two morphisms / : A — > V* 
and g : A — > A* such that g is prolongeable on a letter a. The question is to decide 
whether or not the infinite word f{g u) (a)) — W0W1W2 ■ ■ ■ is ultimately periodic. 
Thanks to [36], one can canonically build an abstract numeration system S = 
(L, £, <) and a deterministic finite automaton with output M. — (Q, q , 2, 8, T, r) 
where r : Q — > T is the output function such that 

Vn > 0, w n = T(S(q ,vep s (n))). 

Such a sequence is said to be an S '-automatic sequence. Notice that (w n ) n >o is 
ultimately periodic if and only if for all b £ T, the S- recognizable set 

X b = {n | w n = b} 

is ultimately periodic. If / and g are such that the associated numeration system 
S satisfies the assumptions of Theorem [52] then one can decide whether or not Xb 
(and therefore (w n ) n >o) is ultimately periodic. 
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